Mechanical and superfluid properties of dislocations in solid ^He 
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Dislocations are shown to be smooth at zero temperature because of the effective Coulomb- 
type interaction between kinks. Crossover to finite temperature rougnehing is suggested to be 
a mechanism responsible for the softening of *He shear modulus recently observed by Day and 
Beamish (Nature, 450, 853 (2007)). We discuss also that strong suppresion of superfuidity along the 
dislocation core by thermal kinks can lead to locking in of the mechanical and superfluid responses. 
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Network of superfluid (SF) dislocations [ij is a 
most likely scenario for the supersolid response 
of solid ^He on rotation 0] at temperatures T < 
Q.2K, at least, for values of the SF fraction Ps{T) < 
0.01. The proof [1] that ideal crystal of ^He can- 
not be a supersolid [1] and the observations of SF 
cores of some dislocations in first principle sim- 
ulations 01 put the model Q on a sohd ground. 
Dislocations are building blocks of most topologi- 
cal defects. Thus, understanding their properties 
is of primary importance. 

Classical and quantum mechanical behavior of 
dislocations was addressed by many investigators 
in the past Quantum roughening of disloca- 
tions has been proposed to be important for in- 
ducing supersolidity Q. Remarkable resemblance 
between shear modulus temperature dependence 
G{T) and Ps{T) was uncovered recently [8|. 

Here we study quantum behavior of a single dis- 
location and its crossover to classical regime. We 
also discuss a possibility of suppression of SF along 
the core Q by geometrical kinks on dislocation 
line. Our findings are that such simple model, 
which ignores collective effects of the dislocation 
network, can naturally explain the key features of 
the experiment Q . Some of our preliminary results 
have been presented in Ref.0]. 

Model. Edge dislocation moving along its glid- 
ing plane {x,y) is modeled as a string character- 
ized by some displacement field j/(x, t) with kinetic 
and tension energies Hg = J dx[{ni/2){y tvY + 
{y^ni/2){y xuY] and subjected to Peierls poten- 
tial C/p(y) = J dxucos(2TTy(x,t)/b) with some am- 
phtude u (see in Ref.jlOj). Here ni is Id mass 
density; Vd - speed of sound along the string; 
b denotes Burgers vector. The Hamiltonian Hg 
alone was successfully used for describing dynam- 



ical properties of many materials at high temper- 



atures [llj. At low T, however, the Peierls po- 
tential is essential. Furthermore, static kinks in- 
teract with each other through effective Coulomb 
potential Vcix) 0, 12 1 (induced by exchanging 
bulk sound phonons with bulk velocity Vb ~ Vd). 
There is also long-range contribution to the effec- 
tive mass [l3|- So, we modify Hs + Up by including 
the long-range interaction with retardation effects 
and obtain Vc from the bulk phonon propagator 
integrated over the directions of momenta perpen- 
dicular to the core. Thus, full action S in imagi- 
nary time is given as 
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in units of typical cell size xq k, 4A of ^He. Unit 
of time is tq = x^jVd so that /? = Tq/T, Tq ~ 
hVd/xi) « bK for Vd ~ 200m/s; and expan- 
sion is performed over Matsubara frequencies uj = 
(27r//3)n, 71 = 0, 1,2, ... and spatial wavevector q^ 
along the core of length = N^xq ; qo ~ 1 
K = Trh/{4niVdb'^); we have chosen the rescaling 
y = 2b6/y/TT, with 9^^g^ standing for the Fourier 
component of 0; v = (Vd/Vb)'^ (we will be using 
V = 1). C describes relative energy of the bulk de- 
formations with respect to the core. The first one 
is usually the largest l^, 13]. Thus, in real mate- 
rials C > 1. In calculating the partition function 
Z = j Dy cx-p{—S/h) we use periodic boundary 
conditions 6'(x-|-i,t) = e{x,t), 0{x,t+(3) = 9{x,t). 
It is important to note that, for solid "^He , the 



Luttinger parameter K in the term ([T]) is close to 
unity, that is, for C = the dislocation is close 
to Berezinskii-Kosterlitz-Thouless (BKT) transi- 
tion (see in Ref.[l3l), so that it is in quantum rough 
state for K > 1/2. Below we will show that arbi- 
trary small C inevitably brings the system into the 
quantum smooth state as L^, ^ oo for any K . 

Following Granato-Liicke change of G can 
be related to a response of a typical dislocation 
segment of the dislocation network. If the whole 
network is viewed as a set of blocks of free seg- 
ments of sizes Lx,Ly,Lz along the correspond- 
ing orthogonal displacement y{x,t) un- 
der the force / = azyb, where Uzy stands for the 
stress tensor [ij], results in a strain of the block 
Uzy ~ b{y) /{LyLz). Thus, (y) cx azy and, accord- 
ingly, Uzy cx; <7zy This leads to 



Given the low w, qx asymptotic of Vc in Eq.Q, the 
space-time asymptotic can be found as Ucaif) ^ 
ln(l + 2C\n.r)/C ~ Inlnr, r = y/WTx^ . That 
is, slower than ~ Inr. Since the entropic con- 
tribution scales as ~ — In?', the free energy of a 
pair of charges UcGif) — Inr ^ — oo becomes un- 
bounded from below. According to the Kostcrlitz- 
Thouless argument this implies deconfinement of 
pairs of charges for arbitrary small C, that is, the 
smooth state without zero-point kinks. 

RG argument. We construct RG flow equations 
in the one- loop approximation for u, K and G by 
considering them as scale dependent. Then, intro- 
ducing the rescaling variable /, we find 
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where Gei stands for the elastic modulus of ideal 
crystal; Ud = h"^ /{LyLz) ~ b^/L^ is the dislocation 
density (in units of b), provided Ly fa Lz ~ L^- At 
large T, where the Peierls potential is irrelevant, 
the system is Gaussian and the high-T modulus 
Goo = G{co,Lx) obeys 
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where g/, = 27r/ L^ for chosen boundary conditions 
and, in general, qi^ cx l/L.^., so that the correc- 
tion to the modulus is only weakly dependent on 
Lx'- (X 1/(1 + 2G\n{Lx/b)). Absence of quantum 
roughening (proved below) implies that the Pcicrls 
potential is always relevant at T = so that the 
linear response is dominated by the lowest term 
in the expansion of cos(...) in Eq.Q. Thus, for 
Go = G{0,Lx) we find 
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For Lxjb » ^jGeib'^/u/{2n). At finite T, the 
dislocation undergoes a crossover from quantum 
smooth, Eq.®; to classically rough state, Eq.([5|). 

Heuristic argument against quantum roughen- 
ing. The system (jll2|3p can be mapped on the 2D 
classical gas model characterized by the interaction 
(in Fourier) between the integer charges which are 
dual variables to kinks Il3l: 
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An elementary analysis shows that, regardless of 
the initial conditions, the flow is always toward 
K{1) 0,u{l) oo, that is, to the gapped 



state. This conclusion is similar to the result [15 1 
for the case of dipole-dipole interactions in Ising 
model which also maps on Sine-Gordon model. We 
note, however, RG equations (|8l9ll0p have only 
qualitative value and, as a comparison with the 
Monte Carlo (MC) simulations show, do not de- 
scribe quantitatively the actual flow. 

MC simulations. Treating the term ([2]) in the 
Villain approximation (see in Ref . [l3| ) . the action 
(|1I2I3P is reformulated as 



5 = i ^ Ud{x.i - Xj)JiJj 
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in terms of the conserving integer currents J; = 
J^"\ V = ±a;,±f, deflned on each bond coming 
out of a given site i of the space-time lattice, with 
zbi, ±i being respective unit vectors. The Fourier 
of Ud{x) is 
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with a = uAxAt and A^,, A( being respective 
space-time sizes of the discretization cell (in units 
a;o,To). Eq.(|4]) in the J- variables becomes: 
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FIG. 1: (Color online) Stiffness K{N^,C) at T = as 
a function of ^ = Cln A^^;. Shown numbers are values 
of C. Chosen parameters are K = 1.6, a = 0.01. Inset: 
family of curves K{Nx, C) for various C. 



FIC. 2: (Color online) Experimental shear modulus 
G(T) from Ref.[|| and its fits by Eq.jni for three dif- 
ferent sets of parameters, A — 0.0888. Error bars are 
smaller than symbols sizes. 



in the limit Lr,. oo. Here 7 = b^naL^GQ — const; 

Wt = J2i jf'' /Nt, Nt = /3/At is the total current 
in the t-direction, so that k is the compressibility 
of the J-current model (jll|) . 

MC simulations of the action (fTl|) have been 
conducted by the Worm Algorithm [l6|. We used 
isotropic lattice Nt = Nx with A^; = 1 and At = 1, 
in progression of sizes and measured func- 
tion of iV^, C for A' = 1.6 (well above the BKT 
critical point K = 0.5 for C = 0). The renormal- 
ized stiffness at q^. = 2'k/Nx for Nj. >> 1 is defined 
as K{Nx, C) = K- (1-|-Vc(9k, 0)). For given param- 
eters and various values of C all the data (shown 
in the inset of FiglT]) can be collapsed on a single 
master curve K{Nx, C) = K{£,) where ^ = ClniVj, 
and K{£,) = K{£^, 1), FiglU which asymptotically 
reaches zero. This concludes our proof that a sin- 
gle dislocation is in the quantum smooth state for 
an arbitrary value of the long-range interactions C. 

Finite T behavior. At T = there is a gap 
A in the spectrum of normal excitations. Thus, 
T-dependence G{T) of G becomes significant only 
at r > A. As T increases further, k becomes fi- 
nite and eventually reaches the asymptotic value 
= k{oo,Lx) defining Eq.®. Specific shape of 
k{T,Lx) has been obtained from the MC simula- 
tions of the action pi|) at various T ~ To/{AtNt) 
for fixed (At varied as 0.5 — 1.5 and it has been 
checked that the result does not depend on the 



time discretization). As it turns out, k can be fit 
by some master curve as well: k/koo = F{T/Ta), 
with F[x) being some universal function; depen- 
dence on the parameters C,K,a,Nx, is hidden in 
Ta. Thus, for any T Eq.([T3D becomes,Figll 

Core superfluidity and shear modulus stiffen- 
ing. One of the most striking results of Ref. [8| is 
the similarity of the temperature dependencies of 
Ps{T) and G{T). Here we demonstrate that such 
feature can naturally arise if the core SF is strongly 
suppressed in the presence of the geometrical kinks 
formed on segments of the dislocation network at 
finite T. A typical T-scale where k becomes fi- 
nite is given by A so that it is natural to expect 
that SF is suppressed at T = Tc ~ A rather than 
at some intrinsic for smooth dislocation, pro- 
vided A << Tc- Phenomenologically, this implies 
the term &' «;(T)|?/jp, where 6' > is some coeffi- 
cient in Landau free energy i/^, with standing 
for an effective 3D SF order parameter. The SF 
part can be written as — 6o(T)|V'P + i3|V'l'*/2, with 
some constant B > Q and function b(){T) < for 
T > Tc and boiT) > for T < Tc. If fc « Tc, 
one can set T = in b(){T) and obtain 

i^L = foo(0)(^-l)|^P + ||^|^ (15) 



where 6' has been ehminated from the condition 
that the actual transition occurs at T^- Thus, using 
p^(0) = |^/;(0)p = &o(0)/B (since k(0) = 0), one 
finds 
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for T < and p^(r)/p^(0) = otherwise. 
We define relative variation of the modulus (fT5|) : 
g{T) = (G(T) - G(f,))/(G(0) - G(f,)) and find 
from Eq.p^ 
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So, 



within 10% of accuracy (since jk{T) < 0.1) 
Eas. lfBTTTD imply p,(r)/p,(0) « ^(T). 

Summarizing, we have shown that a single dis- 
location must be smooth at T = 0. Crossover to 
finite temperatures leads to classical roughening 
and, accordingly, to intrinsic softening of the shear 
modulus. Our model contrasts with the model 
0, where the central role in the softening effect 
is endowed to the '^He impurities boiling off from 
the dislocation cores and, therefore, eliminating 
pinning centers. Within a minimal model, where 
superfluidity along dislocation core is suppressed 
by geometrical kinks, it is possible to understand 
the similarity between SF and dynamical responses 
0, i of solid "He. 

Present study focuses on a simplified model 
where collective effects of the dislocation network 
are ignored and, therefore, plenty of unanswered 
questions remains: What is the role of finite den- 
sity of dislocations, in general, and, specifically, 
can there be a screening of the long-range interac- 
tions for certain type of deformations? Can a ID 
crossover to classically rough state at finite T be- 
come an actual phase transition in the 3D network? 
What is the behavior of slanted dislocations? Our 
model does not address the issue of "^He impuri- 
ties. They provide non-periodic trapping potential 
for dislocations and, therefore, their role in modi- 
fying G(T) requires a separate analysis. 
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